Numerical Analysis and Scientific Computing Preprint Seria ILU preconditioners for non-symmetric saddle point matrices with application to the incompressible Navier-Stokes equations

نویسندگان

  • I. N. Konshin
  • M. A. Olshanskii
  • Yu.V. Vassilevski
  • IGOR N. KONSHIN
  • YURI V. VASSILEVSKI
چکیده

Motivated by the numerical solution of the linearized incompressible Navier–Stokes equations, we study threshold incomplete LU factorizations for non-symmetric saddle point matrices. The resulting preconditioners are used to accelerate the convergence of a Krylov subspace method applied to finite element discretizations of fluid dynamics problems in three space dimensions. The paper presents and examines an extension for non-symmetric matrices of the Tismenetsky–Kaporin incomplete factorization. It is shown that in numerically challenging cases of higher Reynolds number flows one benefits from using this two-parameter modification of a standard threshold ILU preconditioner. The performance of the ILU preconditioners is studied numerically for a wide range of flow and discretization parameters, and the efficiency of the approach is shown if threshold parameters are chosen suitably. The practical utility of the method is further demonstrated for the haemodynamic problem of simulating a blood flow in a right coronary artery of a real patient.

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

ILU Preconditioners for Nonsymmetric Saddle-Point Matrices with Application to the Incompressible Navier-Stokes Equations

Motivated by the numerical solution of the linearized incompressible Navier–Stokes equations, we study threshold incomplete LU factorizations for nonsymmetric saddle-point matrices. The resulting preconditioners are used to accelerate the convergence of a Krylov subspace method applied to finite element discretizations of fluid dynamics problems in three space dimensions. The paper presents and...

متن کامل

Numerical Analysis and Scientific Computing Preprint Seria An algebraic solver for the Oseen problem with application to hemodynamics

The paper studies an iterative solver for algebraic problems arising in numerical simulation of blood flows. Here we focus on a numerical solver for the fluid part of otherwise coupled fluid-structure system of equations which models the hemodynamics in vessels. Application of the finite element method and semiimplicit time discretization leads to the discrete Oseen problem on every time step o...

متن کامل

A comparison of preconditioners for incompressible Navier-Stokes solvers

We consider solution methods for large systems of linear equations that arise from the finite element discretization of the incompressible Navier–Stokes equations. These systems are of the so-called saddle point type, which means that there is a large block of zeros on the main diagonal. To solve these types of systems efficiently, several block preconditioners have been published. These types ...

متن کامل

Efficient Solution of Large-Scale Saddle Point Systems Arising in Riccati-Based Boundary Feedback Stabilization of Incompressible Stokes Flow

We investigate numerical methods for solving large-scale saddle point systems which arise during the feedback control of flow problems. We focus on the instationary Stokes equations that describe instationary, incompressible flows for moderate viscosities. After a mixed finite element discretization we get a differential-algebraic system of differential index two [J. Weickert, Navier-Stokes Equ...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:

دوره   شماره 

صفحات  -

تاریخ انتشار 2015